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Abstract 

The purpose of the present article is to call attention to some realistic quasiparti- 
cle-based description of the quark/gluon matter and its consistent implementation in 
thermodynamics. A simple and transparent representation of the thermodynamical 
consistency conditions is given. This representation allows one to review critically 
and systemize available phenomenological approaches to the deconfinement problem 
with respect to their thermodynamical consistency. A particular attention is paid to 
the development of a method for treating the string screening in the dense matter 
of unbound color charges. The proposed method yields an integrable effective pair 
potential which can be incorporated into the mean-field picture. The results of its 
application are in reasonable agreement with lattice data on the QCD thermodynam- 
ics. 
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1 Introduction 



With the advent of RHIC and LHC, there is a growing need for a deeper understanding 
of various properties of the QCD matter at high temperature and finite density. At the 
moment, we are still far from a satisfactory level in this respect, even for equilibrium 
properties of quark-gluon plasma. Indeed, though such a system can in principle be 
approximated as a gas of quarks and gluons, a fully perturbative calculation with these 
degrees of freedom does in practice not work well at any reasonable temperature since 
the perturbative series are badly converged due to infrared-sensitive contributions. On 
the other hand, the QCD lattice calculation, the only systematic fully nonperturbative 
method available, is restricted in the presence of light dynamical quarks, and even more 
so in the presence of a finite baryon density (see [|l| where the current state of art is 
summarized). Therefore, various phenomenological, QCD- motivated models are called 
up for describing the thermodynamics of highly excited nuclear matter and its Equation 
of State (EoS). 

General arguments from QCD and lattice data tell that a kind of string is developed 
between quarks and antiquarks at large distance and it is natural to identify such qq 
system with conventional mesons. Treating quark and gluon propagation in the confining 
QCD vacuum within non-Abelian SU(3) gauge theory, the string dynamics was success- 
fully applied to conventional mesons, hybrids, glueballs and gluelamps. However, if such 
a string is surrounded by unbound quarks and gluons, the qq system can be excited not 
only in color-singlet states, but also in color-octet states or even dissociate into constituent 
elements. The latter will signal, in general, on the deconfinement phase transition. These 
phase transformations are intimately related to the change of string properties (in partic- 
ular, color charges of quarks are screened in quark-gluon environment): string behavior 
becomes medium dependent. 

By now, there is a number of simplified models for describing static hadron prop- 
erties as well as a highly excited, deconfined state of quark matter, the quark-gluon 
plasma (QGP). A common feature of these models is that they all are based on a quasi- 
particle picture, considering isolated particle-like degrees of freedom and assuming that 
these quasiparticles are moving in a background mean field. Two- and many-particle 
correlations are included in the mean-field contributions and in the modification of the 
one-particle spectra. Well-known examples are the original bag model and its later ver- 
sions 0, H, ^, phenomenological approaches with temperature-dependent bag constant 
0, ^, string-motivated density-dependent corrections to an ideal (massive or massless) 
quark matter equation of state, the very consideration of hadrons as composite objects 
in QCD, excluded volume corrections 0, |^, and finally mixed phase and chemical 
mixture |T0| , |Tl| models dealing with the transition between quark matter and hadron 
matter in a phenomenological way. 

The present paper concerns the quasiparticle description of the QCD thermodynam- 
ics with the particular emphasis on the mean-field treatment of in-medium strings. The 
paper is organized as follows. In Sec. II we consider the thermodynamical consistency 
of the quasiparticle description in general. Any phenomenological approach, involving 
a quasiparticle interaction, usually operates with a Hamiltonian which may depend on 
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thermodynamical characteristics of the surrounding matter, hke the temperature T and 
density n. As is known for a long time (see, for details, [0, |13|), there exist certain 
restrictions to the dependence of such a Hamiltonian on the thermodynamic variables. In 
this section a transparent and useful representation of these restrictions is derived [see 
Eqs.([21|) and (p2D] which directly involves the quasiparticle spectra. The obtained repre- 
sentation of thermodynamical consistency allows us to get an instructive relation between 
a number of modern approaches, dealing with the deconfinement problem in the frame- 
work of a quasiparticle picture, as exemplified in the end of Sec. II. The QCD-motivated 
interactions, in particular string-like interactions [|14|, are under investigation in Sec. III. 



A comprehensive model of string formation in the dense matter of unbound color charges 
is developed, which supports the choice of the mean field proportional to an inverse power 
of the color-charge density, as proposed in the papers 0, 0, 0, . This quasiparticle 
scheme is applied in Sec. IV for thermodynamics of the deconfined QCD phase. The case 
when the system is in thermal equilibrium but not in chemical one is also considered. The 
results are summarized in the concluding Sec.V. 



2 Quasiparticle Hamiltonian 



In order to obtain an effective quasiparticle description of a medium made of unbound 
colour charges, one should operate with screened long-range potentials. A natural way to 
introduce the screening in quark matter is based on using the probability density P(r) to 
form a string of the length r. It is worth noting that this investigating scheme has much 
in common with another one which deals with the probability density that the nearest 
neighbour occurs at a distance r [|T6|. Both approaches involve thermodynamical variables, 
which lead to a screened pair potential depending on thermodynamical quantities. Due to 
such effects the quasiparticle Hamiltonian becomes density and temperature dependent, 
which in turn leads to a modification of a thermodynamical potential (e.g. the Gibbs free 
energy). Eventually a nonideal EoS emerges. 



2.1 General structure of Hamiltonian 

We start with the general quasiparticle Hamiltonian 

H = Y.Y.^ki{T,n) aLak. + l^$(T,n). (1) 

i k 

Here a\^^ and Oki are the usual creation and annihilation operators for quasiparticles of the 
i-th sort with momentum k. They also may depend on other internal degrees of freedom 
like spin, color, isospin, etc. The volume V is constant and large enough (infinite in the 
thermodynamical limit). So in this limit, the summation over quantum states labeled by 
k can be replaced by a phase space integral (including the internal degeneracy factor di): 
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In general, the one-particle energy eki(T, n) and the background field contribution $(T, n) 
to the energy density, both depend on the temperature T and the set of particle densities 
n = {rii, n2, . . .}. Notice that the quantity V^(T,n) is nothing else but the energy of 
the quasiparticle vacuum. Generally speaking, it differs from the vacuum of primordial 
particles, which leads to the c- number term appearing in Eq. (^. In the present paper 
we consider the situation when, similarly to the case of the Hartree-Fock quasiparticles, 
the expectation value of the quasiparticle number operator iVj = J2k '^L'^ki equals to the 
number of primordial particles Ni = riiV. This implies that we deal with the picture of 
quasiparticles interacting and, thus, correlating with each other. In turn, the expecta- 
tion value of the Hamiltonian has to be equal to the mean energy of the system under 
consideration. This leads to the following relations: 

<H> = E = Vj2diJ j^ekiiyM + V^, (3) 

with the occupation numbers i/^j = (aj^-Oki). There is another way of calculating the mean 
energy E = Ve{T,n) and mean multiplicity = Vui which proceeds from a thermo- 
dynamical potential rather than from Eqs. (|^) and (^. For density- and temperature- 
dependent Hamiltonians these different ways may lead to different results (see, for ex- 



ample [y, 12|). Hence, in what concerns the dependence on n and T, the quasiparticle 
Hamiltonian will have a correct structure only if the thermodynamical consistency re- 
quirements and @) are satisfied when starting with either the Hamiltonian or the 
thermodynamical potential. 

2.2 Chemical potentials 

Using temperature T and number densities rii = Ni/V as basic descriptive variables, the 
thermodynamical behavior and the appropriate EoS can be derived from the correspond- 
ing thermodynamical potential, the free energy F{V, N,T): 

F^Vf = TY^^,d^l -0^ In (l - ^^ e-^-) +Vj2fi,n, + V^. (5) 

Here = ±1 is determined by the quasiparticle statistics, /ij stands for the chemical 
potential of the quasiparticles of the i-th sort, and Zki is defined by 

Zki = • (6) 

The simplest way of calculating the free energy implies the use of the grand canonical 
ensemble when particle numbers are known only in average and chemical potentials, /Zj, 
are introduced instead of Ni as descriptive thermodynamical variables. To calculate the 
partition function in this case, the quasiparticle Hamiltonian (|I]) should be modified to 

H' = H-Y.^i^N^, (7) 
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where H is defined by Eq. (|T]). We recall that the physical meaning of /Xj is the energy 
loss by removing a quasiparticle of the i-th species while the total entropy and volume 
of the system are kept constant. This chemical potential a priori has nothing to do 
with the fact whether this particle really carries a conserved charge or not. However, 
there are Lagrange multipliers associated to the conservation laws of such charges like 
the baryon number, strangeness or electric charge. In order to elucidate the difference 
between chemical potentials in general and those associated to the conserved charges, let 
us consider a particle mixture of many sorts whose abundance is known only in average. 
The mixture components (not necessarily all of them) carry some conserved charges. The 
conservation of these charges is controlled by the appropriate chemical potential. We 
denote such a charge of type h carried by a particle belonging to the i-th component of 
the mixture as qu- Then for conserved quantities we have 

i 

Usually there are more components than the number of conserved charges. In particular, 
it is the case for quark-gluon matter, to which we pay special interest in the present 
paper. Besides the case of a one-component system, the set of Eq. (|^) is insufficient 
for calculating all mean numbers for the mixture components. Therefore, we need some 
additional requirements which would allow us to determine the particle numbers Ni by 
making use of Eq. (|]). The chemical equilibrium is usually assumed and then these 
additional requirements are given by 

/^r = n (9) 

b 

where stands for the chemical potential associated to the charge sort h. In the general 
case, when the system is out of chemical equilibrium and the component concentrations 
become time-dependent, the chemical potential can be split into two parts, 

A^i = J2^ibfJ'b + (ii , (10) 
b 

and Eq. (0) is reduced to 

H' = H-Y.^lbQb-Y.J^^^^■ (11) 
6 i 

The quantities /Ij describe the deviation from chemical equilibrium in the thermally equi- 
librated system. They are exactly zero in the chemical equilibrium limit, resulting in the 
familiar relation 

K,=H-j:fibQb. (12) 

b 

Below we shall investigate consistency of the quasiparticle picture in thermodynamical 
treatment including the possibility of deviations from chemical equilibrium in a mixture. 
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2.3 Thermodynamical consistency 

As mentioned above, any approach starting with a thermodynamical potential appears 
to be thermodynamically consistent. In other words, if all the quantities of interest can 
be calculated only through the derivatives of this thermodynamical potential, one is pre- 
vented from encountering thermodynamical inconsistency. Problems arise, however, when 
calculation can proceed not only from the constructed thermodynamical potential but also 
from a more fundamental level, some quasiparticle Hamiltonian at a given temperature 
and/or density. In the last case the result may depend on the calculation method un- 
less the quantities eki(T,n) and $(T, n) obey relations derived in accordance with the 
consistency requirements (|^) and (^). These relations are below called as conditions of 
thermodynamical consistency. 

To elaborate on these conditions, let us consider a system with the Hamiltonian defined 
by Eqs. (|l|) and (|^). If the abundance of all mixture components is known only in average, 
the proper thermodynamical potential has the form 

n = E-TS-J2fi,Ni, (13) 

i 

with 5* being the total entropy. The grand-canonical partition function corresponding to 
H' can readily be calculated 

Z = tr(e-^'/T)^ g-w/T^g-n/r_ ^^4^ 

\k,i nki J 

Recall that in the Fermi-g the integer quantity n^i equals either or 1, while in 



the Bose one it runs from to 00. Hence, from Eq. (|T4D it follows that 



^=ri:<;.{./^ln(l-&e-".)+4. (15) 

Note that Eq. (|15D covers Eq. (1) with the definition F = E — T S . It is well-known that 
for the Hamiltonian H' the average quasiparticle occupation number u^i is 

^ki = {nki) = ^ , (16) 

which together with Eq. ([T5|) leads to 

r d^k 

- = - Tj2dSj In (1 + z/fc.) . (17) 

The total differential of the thermodynamical potential Q is given by 
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Since in the grand-canonical ensemble the quasiparticle densities Ui and average occupa- 
tion numbers z/^j are functions of the temperature T and chemical potentials all the 
differentials can be expanded in terms of dT, dV and d^i. In particular, we have 



duki 
d^ki 
d^ 



-m (1 + ii m) d 

deki , deki 

dfjLj + di , 



(-ki 



T 



(9$ 



dT 

a/^i + -gf dT. 



Inserting this relations into Eq. (piq ) and comparing then the derived result with the 
general formula 

dn = -p dV - S dT dfii , 



one can arrive at 



di 



d^k 



(27r)3 
d^k 



^ki + $ /iiCi - TCt ■ 



where 



(27r)3 
df 



^ki ~ Ci 



(19) 
(20) 



d'^k dekj 



9$ 



(27r)3 dT 
d^k dek- 



(27r)3 dfii 



^kj-i 



Equations (|1^) and (pOf ) should be compared to the consistency requirements given by 
Eqs. (|^) and (^). It leads to the conditions Ct = and Ci = 0. 

This result can be rewritten in a more manageable form which allows for relating our 
specific case based on Eq. (|l|) to a more general one. To elucidate this connection let us 
consider the derivative matrix 

drii 



M, 



Its elements can, in principle, have arbitrary values, and we expect that the determinant 
of Mij is not zero. This is indeed the case since the derivative matrix is given by the 
following implicit relation: 



T 



di 



d^k f dei 

Vkj (1 + Z^fci) Oij 



(27r 



leading to the matrix equation 
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with 



1 f (Pk 
Ajk = 6jk — dj J Ukj (1 + Ukj), 

1 f d?k dek 



^ij = f ^3 J (1 + l^kj) 



(27r)3 ^ driki' 

This matrix equation has a formal solution M = (1 + 5)^^ A. While i? may have zero 
eigenvalues, A does not, so the determinant detM = detA/det{l + B) cannot vanish. 
Using this information, the equalities Ct = and Ci = can be rewritten as follows: 

9$ f d^k deki , . 

a^ + i:^./p^a^-^. = 0. (22) 

Equations (^T]) and (^21) represent a particular case of the more general relations f^, |T2l 

(^)^O. (^?g£).o. ,23) 
These conditions of the quasiparticle consistency are reduced to Eqs. ( ^TD and ( ^21) when 



the (temperature- and density- dependent) effective Hamiltonian -ffe// has the quasiparti- 
cle form (p. Equations ( ^31) have first been derived in [1^ (see also Ref. 0) under the 
requirement that all the statistical ensembles of the system governed by the density- and 
temperature-dependent Hamiltonian -ffe// yield the same thermodynamics in the infinite 
volume limit. Thus, to avoid thermodynamical inconsistency, the constraints ( P5D should 
be satisfied while constructing the effective quasiparticle Hamiltonians. 

Let us emphasize that thermodynamical consistency is not sufficient by itself when 
thermodynamics is constructed starting from the level of a thermodynamical potential 
but the Hamiltonian structure is ignored. In this case nonphysical expressions can be 
involved even if there is no problem with thermodynamical consistency and all the ther- 
modynamical quantities are derived by differentiating a thermodynamical potential. We 
suggest that relations (|21| ) and ( ^2]) should also be employed in the situation like that to 
avoid unreasonable expressions for quasiparticle spectra which can be met in the literature. 



For instance, see the papers on the compressible bag model |18|. It is thermodynamically 
consistent but the quasiparticle spectra used there have nothing to do with (^) and (^2]). 
Another example concerns the approach of Ref. that has no problem with thermo- 
dynamics, too. However, the shift of the chemical potential proposed in that article is 
equivalent to the introduction of a temperature-dependent vector-type mean field. It is 
shown below from Eqs. (^) and (^) (see the next section. Example 1) that the mean 
field like this can not depend on the temperature explicitly. By the way, it is quite pos- 
sible that the nonphysical feature of the quasiparticle spectra used in the papers is 
an actual reason for causality violation when the sound velocity is getting larger than the 
velocity of light (for more details, see Ref. 
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Note that when temperature and density dependence of the effective Hamiltonian 



is mediated only by some thermodynamical quantity A, the consistency conditions (|23|) 
are equivalent to (dHeff/dA) = 0. This relation comes from the well-known expression 
5F = {6Heff) ||2l|, where 5F and SHgff stand for infinitesimal changes of the free energy 
and Hamiltonian, respectively. This expression can easily be derived by analogy with 
the familiar Hellmann-Feynman theorem and, taken in conjunction with the extremum 
condition for the free energy with respect to the parameter A, leads to (dHeff/dA) = 0. 

2.4 Quasiparticle spectra 



Conditions of thermodynamical consistency (0) and (22) result in certain physical re- 
strictions to the mean-field potential depending on the structure of quasiparticle spectra 
without coming into any detail of interaction between constituents. We shall demonstrate 
that in a few cases used in phenomenological treatments. 

Example 1. Let the energy of a quasiparticle of the i-th sort moving with the 3- 



momentum k be approximated as (for instance, see Refs. [T^, [T^, |T^, ^ ^ 



ek^{T,n)=Ui{k) + Ui{T,n), (24) 

where uJi{k) stands for the energy of the free particles of the i-th sort and, as above, 
n denotes the set of particle number densities {ni,n2, . . .}. Generally, the mean-field 
potential Ui{T,n) is a function of temperature and particle densities. The free particle 



energy can be given in either the relativistic form, uji{k) = yA;^ + rriQ^, or in the non- 
relativistic one, uJi{k) = k'^/2moi. It depends only on the momentum k and bare particle 



mass moj. Then, the conditions of quasiparticle consistency (^) and ( p2D are reduced to 
the following equations: 

— + > n, — ^ = 0, h> n, — ^ = 0. (25) 

dT ^ ^ dT ' diii ^ ^ dui ^ ' 

As it is seen, ignoring the background field contribution $(T, n), like, for example, in the 
papers |T7| , results in a loss of thermodynamical consistency. It is important to note 
that Eqs. (^) are only then compatible, if the mean field Ui does not depend explicitly 
on the temperature. Indeed, the first equation in (|25|) has an integral of the form 

E ^jUj + $ = 
j 

where ip = (p{n) is an arbitrary function of the quasiparticle densities. By differentiating 
this expression with respect to n^, we get 

— -f/+Vn ^£i_^^ 



Taken in conjunction with the second equation in (|25D , the obtained relation is reduced 
to dip /drii = Ui. It follows then that Ui and $ are temperature-independent functions. In 
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other words, when quasiparticle spectra are defined by Eq. (p^, the thermodynamically 
consistent mean-field potential may depend only on particle densities: f/j = f/j(n), $ = 
Note that this important point is missed in some papers |25|, where the excluded 
volume effects are treated by means of the mean-field approximation. 
Equations (|25| ) lead to one more interesting result 



m 

drij 



dm 



(26) 



This crossing relation, first presented in |^, follows from the second equality in Eqs. (^Sf ). 
To derive Eq. (p6|), one should differentiate the second equation in (^) with respect to 
rii. Then, by interchanging the indices / and i and comparing the obtained expression 
with the previous one, we arrive at Eq. (EB). By doing so one should keep in mind that 



52$ 



52$ 



drii drii drii drii 



drii drii drii drii 



This is valid, provided the second derivatives of Uj and $ are continuous functions of T 
and n, which is usually the case. Note that Eq. (^) is very useful when dealing with the 
mean fields f/, for many-component systems. For example, see the investigation of quark- 
hadron interactions in 0. The crossing relation (|26| ) should be kept under constructing 
the mean fields acting on quasiparticles of different species in a many-component system. 
Otherwise, the thermodynamical consistency can be lost. As an example, one can point 
out Ref. |2^, where mean fields were chosen as Ui oc ritot {Tn/niQiY with 5 = 1 or 2. Here 
m denotes the nucleon mass and Utot = H^j- 

Example 2. Another popular form of quasiparticle spectra eki{T,n) is given as 



efc,(T, n) = ^k^ + mKT,n), (27) 
i.e. an effective quasiparticle mass rrii = mi{T,n) is introduced in a way similar to scalar 



mass in the relativistic mean- field theory [^]. In this case Eqs. ( pT]) and (^) can be 
rewritten as follows: 



dT 



n 



drrij 
'df 



0, 



drij , drij 



0, 



(28) 



where the Lorentz-scalar quasiparticle density n[^^ is defined by 

(s) , f d^k rrij 



(27r) 



(29) 



For the sake of simplicity, let us limit ourselves to the case of one sort of quasiparticles. 
Differentiating the first equation in ( P5| ) with respect to n and the second one with respect 
to T, we arrive at 

Qfii^) Qfyi dn^'^^ dm 



10 



provided the mixed second derivatives of $ and m are equal to each other. As follows 
from Eq. (^), the gradients of functions n^'^\T,n) and m{T,n) are parallel vectors in 
the {T,n) plane. Hence, m is left constant along any line where n^'^^ is constant. Since 
n^^\T,n) and ■m{T,n) are different iable functions of T and n, the {T,n) plane is densely 
covered by lines of constant n^'^\ Therefore, in a thermodynamically consistent model 
with one sort of quasiparticles, whose energy is defined by p7|), the effective quasiparticle 
mass depends on temperature and quasiparticle density only through the scalar density. 
The Walecka model [26] without vector field (for zero baryon density) is a particular case 
of the considered variant. Another example can be found in the paper of Boal, Schachter 



and Woloshin in Ref. |jT5|, where interactions in the quark-gluon plasma are described 
by introducing the effective masses of quarks and gluons depending on sum of the color- 
charge densities. As it has been proven above, this version is inconsistent. One can expect 



that the case of many quasiparticle species is described by equations similar to Eq. 
i.e. mi is a function of the set of n\ 
the scalar densities too and satisfies 



i.e. mi is a function of the set of nf\ As a consequence $ depends on T and n through 



$ = - y Y.^f'^^r (31) 

j 

Deconfinement models dealing with the temperature- and density-dependent masses of 
quarks and gluons ^ are also related to the Example 2. Here it is often assumed 

n'^^^ oc m^ and $ oc m^. A purely temperature-dependent bag constant, $(T), without 
mass modifications on the other hand is inconsistent. The same is related to the situation 
when temperature-dependent masses without the background term are used [p9|. To go 



in more detail, see also the papers |]T3| 



Example 3. If the mean field in Eq. ( plf) is scaled with some coupling constant, Ui = 
Qi U (see, for example, Refs. [y, |T^, ^), then conditions of thermodynamical consistency 
take the form 

9$ dU dU 

with 

P = Y.9j ny 

Using the procedure similar to that described in the previous example, one can be con- 
vinced that the density dependence of the mean field U is mediated by p{n) only. If 
Qi is proportional to the baryon number hi of the quasiparticle i: gi = gbi, then we get 
p = gJ2i biUi = gub, where Ub is the total baryon density. Similar situation is realized 
when the quasiparticle interaction is mediated by a vector field. 



Now let us return to Ref. p3[ mentioned at the end of the Example 1. Eqs. (|3^) 
suggest how one can correct the mean field Ui = Utot (jn/mQiY, used in this paper, in such 
a way that to keep the relation Ui oc (jn/moiY. It turned out that the unique solution 
is given by Ui = ip{p) (jn/mQiY, where instead of Ufot we use p = Y.ij^/rnojYnj and an 
arbitrary function ip{p) which can be chosen as ip{p) = p. 

Sometimes it is convenient to subdivide the full set of coupling constants into the two 
groups: gl""^ < (corresponding to the attractive interaction) and gi^'' > (related to the 
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repulsive interaction, see, for example, Ref. |2^). In this case Ui = gf^Ur — \gi^\Ua-, and 
we can expect that the repulsive-interaction component Ur is a function of pr = J^idi^^n-i, 
while the attractive one Ua depends on quasiparticle densities through pa = Y^iQi^^^i. 

Example 4. By analogy with the approximation Ui = QiU considered in the Example 3, 
the effective quasiparticle mass of Example 2 can also be scaled as = gi M{T,n). In 
particular, M{T, n) can be the constituent quark mass, whereas gi is the number of quarks 
in the baryon cluster of the i-th sort. Generally, a cluster of g^ constituents may consist 
of quarks, antiquarks and gluons {i.e. forming baryons, mesons, hybrids, glueballs), as 
well. By doing so, we get 

with 
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The thermodynamics of such a system depends on the descriptive variables T and n 
through the quantity p^^\T,n). 



3 Mean-field treatment of string interactions 

Until now we have discussed general restrictions to the phenomenological Hamiltonian 
due to the thermodynamical consistency. It can be further elaborated by specifying the 
interaction between generic constituents. The quasiparticle properties ought to be derived 
from this underlying interaction. In our phenomenological treatment we consider a par- 
ticular case of strong pair interaction mediated by strings stretched between color charges. 
The main difficulty here is that such a system is plagued with long-range interaction, and 
only in-medium screening renders the problem treatable, even in the weak-correlation 
approximation. 

Strings are particular QCD field-constructions involved in the interaction between two 
color charges in vacuum. We consider here how this interaction will behave in a medium 
consisting of point-like color charges. Let some reference color charge creating a string be 
placed at the origin of coordinates and other color charges be distributed around with the 
density n{t) where ^ is the distance from the reference charge. Physically, we can expect 
that the string formation is characterized by a probability depending on its length, but 
not on its formation history. Let P{i)d£ be the probability for a string to have a length 
between i and i + di. The quantity P{i) can be represented as a product of two factors: 

P{£) =(^1- dx P{x)j w{i) . (34) 

The first factor in Eq. is the probability for the string to have the length not less 
than i. The second factor w{£) is related to the conditional probability, w{i)di, meaning 
that a string is formed between i and i + di (provided it has already reached the length 
i). 
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For gradually growing strings the quantity w{i) can be obtained invoking the argu- 
ments similar to those used in calculation of the mean-free path of a particle moving 
through a medium in a given direction. In this scenario we assume that a string is caught 
by any color charge within a cylinder of the radius a and with axis along the considered 
string direction. This leads to 

w{i) = vra^ n{i). (35) 

The factor vra^ is interpreted as string cross section accounting also for a lack of string 
dynamics. The effective radius, a, may depend on the medium. 

In another scenario the strings are assumed to wildly fluctuate in direction. The gross 
factor w{i) is rather well approximated in this case by the relation 

w{i) = Att f n{e). (36) 

Eq. (|36| ) having an additional £^ factor can be derived by analogy to Eq. (^) if the whole 
area of the spherical shell, Airi'^di, at the distance i is taken into account. Here all possible 
partner charges, located at the distance i, potentially participate in the screening. Note 
that in both cases n{i) stands for the number density of potential partners, on which a 
string of length i can be closed. 

Usually it turns out to be more convenient to deal with the integro-differential equa- 
tions rather than with the integral ones. After differentiating Eq. (|34D, we arrive at 



^ = -wii) Pii) + ^ [l - / P(x) dx ) . (37) 



After substituting here the integral definition (|3^) we get 

§ = -^Q P(i) + § ^ . (38) 

at at w 

This ordinary differential equation is separable and has the following explicit solution 

P{£) =bw{£)e-^o^^^^ '^^ (39) 
The integration constant b is determined by normalization of the probability density: 

/ P{x) dx = 1. (40) 
Jo 

Let us consider now Eqs. (|35D, (^) and (|39D in more detail. Neglecting the spatial charge 
correlations {i.e. taking n{i) = n = const), Eqs. (^) and ( p9|) (belonging to the first, 
straight string scenario) results in an exponential screening 

P{i) = vra^ne"™'" (41) 

/^From here, the probability for a string to be shorter than rj is given by 



Q{r]) = r P{x) dx = l- e"™^" (42) 
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and the average string length is 



(£) = / xP{x)dx=^^. (43) 



Noting that the probabihty distribution (|41| ) can be represented in the form 



P(/)oce-^W/^(<^>) = e-^/<^> , (44) 

one sees that short strings are energetically favored. This looks as a natural conclusion, 
but nevertheless it is not trivial, because "energy" arguments have not been involved 
explicitly in reasoning. 

Within the second scenario of the string screening (i.e. strings are wildly fluctuating, 
see Eq. (p6|)), the corresponding probability density becomes 

P(£) =47rn£V^""^'/3_ ^45) 

This expression covers the result of the papers [|l^ where the probability of string forma- 
tion has been calculated in the nearest-neighbor approximation. The probability of the 
string length to be less than rj has now the form 

Q{r]) = r P{x) dx = l- e-t"'^''', (46) 
Jo 

and the average string length is given by 

= r(-) (47) 



'3 \Ann 

with r(. . .) being Euler's Gamma function. By analogy with the representation (53) 
Eq. (^51) can be rewritten as 



P(£)oc(£)^expl-(r(^)— 1 1, (48) 




which also agrees with the argument that short strings are favorable. 

As noted above the straight-string scenario based on Eq. (|35D involves the effective 
string radius, a, which may depend on thermodynamic variables. Indeed, the string 
survives only in the case when the characteristic length a does not exceed the mean 
distance between neighboring color charges, ro = (47rn/3)^^/'^. Therefore, operating with 
Eq. ( p5D and Eqs. (^lD-(PD, we should employ a < tq. Now, estimating a ^ ctq (c < 1 is 



some constant), we obtain from Eq. (^) 

4 1 ^ ^ 

(49) 

in accordance with Eq. (^71) under the choice c = 2/y'3r(l/3) ~ 0.7 < 1. Summarizing, 
the thermodynamics of string interactions is ruled by the average length of in-medium 



3c3; (7^)1/3 
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strings. Thus, in spite of differences in P{i), both considered scenarios of the string 
screening lead to quahtatively similar thermodynamical pictures. It is interesting that the 
density-dependent interpretation of a in Eq. (IT) leads to a density- dependent screening: 



with the screening mass Mgcr = (3c^/4)^/'^(7™)^/^. At sufficiently high temperatures the 
QCD thermodynamics approaches the Stefan-Boltzmann regime where n ~ T^. This 
yields Mg^r oc T, what is nicely consistent with perturbative QCD. 

The modification of the energy density due to presence of in-medium strings can be 
constructed as 

Ae = na {£) (50) 

with the string tension a. In this case the color constituents of the system are affected by 
the following mean field: 

U = a {i) = A n-\ (51) 

where the constants A and 7 carry information about the sort of color charges and char- 
acter of the in-medium string screening. Note that the results of lattice simulations for 
SU{?>) symmetry can be approximated by 7 ~ 2/3 |T^. This is in qualitative agree- 
ment with our rough estimate 7 ~ 1/3 neglecting the spatial correlations of color charges. 
Indeed, in respect to the thermodynamical character of the EoS, the only fact is decisive, 
that 7 lies between and 1. 

Concluding this section, we sketch how one should, in principle, deal with the case 
when spatial correlations of color charges are taken into account. Let us consider a 
reference color charge placed at the origin. The important point is that the ratio 

M = ^ (52) 

n 



is nothing else but the radial distribution function |3^ which determines the pair particle 
correlations in the uniform system of color charges with the density n given by 

n = lim n{i). 

£^00 

As follows from Eq. (^2|), Eqs. (41) and ( ^5l) operate with g{l) = 1, which corresponds to 



neglecting the spatial correlations of color charges. To go beyond this simplification, one 
should replace n{i) by ng{i) in Eq. (|39D . In particular, using Eq. (^5|), one can derive the 
following equation: 

P{i) = bna^ngii) e-™' "/o ^'(")'^" , (53) 

the average length of the in-medium strings being dependent on the charge-charge spatial 
correlations. Note that for i —>■ 00 the quantity g{i) tends to 1, and we arrive at the 
exponential decay of P{i) with the screening factor which is the same as in Eq. (^T]), 
but with the different normalization constant b. In principle, this difference can lead 
to another estimate of 7 being closer to the lattice result mentioned above. Thus, the 
probability density P{i) should not be identified with the radial distribution function as 



it has been done in the paper |T6 |. 
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4 Application to QCD thermodynamics 



The developed technique allows one to construct a thermodynamical potential in a self- 
consistent way starting from the microscopic level. The Hamiltonian structure of a par- 
ticular phenomenological model is defined by the physical assumptions used. Some simple 
models for the QCD thermodynamics are considered below in order to illustrate conve- 
nience and power of the conditions of thermodynamical consistency (0), ( p2D and to show 
the rationality of our treatment of in-medium strings. 

4.1 Ideal gas in a bag-like model 

Let us consider the ideal gas of particles whose one-particle spectrum is independent of 
medium parameters, 

eki{n,T) = uji{k) . (54) 
Then Eqs. ( PT| ) and (^2p have the form 

leading to a constant background energy $ = 5. It is frequently associated with the bag 
constant. In this relatively simple case the chemical potentials are nevertheless determined 
by the set of the implicit equations 

r d^k 1 



However, in the classical approximation we have Uki « 1 and Eq. ( ^6]) is reduced to the 
expression 

= / e-(-(^)--»)/^ = x.(T) e^»/^. 
Here the second equality defines Xi{T) which relates to the chemical potential as 

^ T .„ (^) . ,57) 

Thus, for the internal energy density we obtain 

£ = 3Tj2ni + B, (58) 

i 

whereas the pressure is given by 

p = TY,ni-B. (59) 

i 

These equations constitute the classical approximation to the familiar MIT bag model 0, 
|3l| being a popular approach of investigating the thermodynamics of the quark-gluon 
plasma. 
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4.2 Temperature- dependent scalar mean field 

The temperature-dependent scalar mean field accounts for a temperature-dependent mass. 
In this situation for the one-component system the quasiparticle spectrum is given by (cf. 



euiT,n) = ^k^ + m\T) , (60) 

the corresponding scalar density (cf. Eq. (|29D) depending only on temperature. The gluon 
and quark plasma (at zero baryon density) is of particular interest where approximately 

m^{T) =ml + g'^T^ . (61) 

Neglecting the derivative of the slowly changing temperature function g{T), from Eq.(|HD|) 



we obtain 

de, m{T) g'T 

= m [T) = • (62) 

Hence, the Ct = constraint leads to 

One obtains for the equilibrium state (yU = 0): 

e = ICT^ + B, p= \lCT^ - B, (64) 
where /C is given by the integral 

lC = d(p-, J_ + (65) 

J (27r)3 gV^^ _^^x' + g' ^ ' 

The (7 = case provides the original MIT-bag EoS. 
4.3 Density-dependent mean field 

Our next example deals with a system of quasiparticles of a single sort with the density- 
dependent spectrum justified in Sec. 3: 

efe(T,n) = cj(fc) + A n"^, (66) 

where A and 7 have been discussed above. In this case Eqs. ( pi] ) and (^) give rise to the 
following form of the background energy density: 

$(n) = A n^-^ . (67) 

1-7 

Hence, in the Boltzmann approximation we find 

^^ = T \n—- + An-\ (68) 
xiT) 
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whereas for the internal energy and pressure one can derive 



"""^ ^ '^'^ coik) e-'^(^)/^ + n'-\ (69) 



xiT)J (27r)3 ^ ^ 1-7 

p = TiT — A n^-^. (70) 

1-7 

As follows from our consideration in Sec. 3, < 7 < 1. In this case the chemical potential 
grows with decreasing density {e.g. due to string pulling) and therefore such systems reveal 
a strong tendency to form clusters. Free sources of strings, or long strings respectively, 
will eventually be purged out of the system. 

4.4 Massless gluons with string interaction 



Let us consider in more detail the case given by Eq. ( |66D for massless SU{3) gluons. The 
choice 7 = 1/3 and A = (2/3)ct satisfies Eqs. (gD and (|lD at a ^ tq = (3/47™)!/^. The 
relation between density and non- equilibrium chemical potential (see Eq. (p!OD) is now 
given by 

_ , [ d^k 1 

7 (2^ e^k+lan-^/^-,)/T _^ ■ 

In the Boltzmann approximation, which is quite appropriate to the system of interest, 
this equation yields 

where x{T) is proportional to T^: 

1/3 

We seek for the solution of the chemically equilibrium state defined by = 0, 

rie, = xiT) e-i'^"- ''/^ . (74) 
This can be transformed into a simple transcendental equation by denoting 



X(T) = (AT)^ A=(^) . (73) 



^ = T </^ ^=Iy^- (75) 



We get 

\z = (f^ . (76) 

This equation has no real solution above the value of a corresponding to a temperature 
^chem- At this chemical critical temperature the l.h.s. and r.h.s. of Eq. (|76D and their 
derivatives should be equal to each other. So, the last condition gives: 

\ = a (f^" . (77) 
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Comparing with Eq. (|76[), we obtain 

= 1/cT, X = ae , (78) 

where e = e^. Finally re-expressing the temperature, we arrive at the result that below 
the value 

Tchem = ^ (79) 

there is no equilibrium solution for the string-like EoS. Assuming 5'[/(3) symmetry for 
massless gluons with = 16 degrees of freedom, we have A ~ 1.175 and arrive at the 
estimate 

Tchem ~ 0.718 , (80) 

which for a = 0.22 GeV^ gives Tchcm ~ 337 MeV. It is noteworthy that a similar relation 
between the color deconfinement temperature and string tension, = (0.6 — 0.65) = 
280 — 305 MeV, has been obtained in the Monte Carlo simulation of the lattice SU{3) 
quenched QCD |T], 

In Fig. |1| the quantities e/T^ and 3p/T^ are plotted as functions of the temperature for 
the system of gluons with the spectrum Eq. (^) beyond the Boltzmann approximation. 
As it is seen, Eq. (^0]) indeed provides a good estimation for the limiting temperature 
Tchem, which is now 303 MeV. The deviation of the pressure and energy curves from 
the ideal-gas value reflects essential attraction even at ~ 2Tc. It is now interesting to 
clarify to what extent our treatment of the in-medium string interactions agrees with the 
nonperturbative lattice QCD. This can be understood with the help of the special quantity 
{e — 3p)/T'^ that is often called the interaction measure. This quantity is directly related 
to remnant interactions that survive in the high-temperature QCD phase because for the 
ideal massless quarks and gluons e = 3p, as mentioned above. In Fig.|^ the interaction 
measure for the SU (3) gluon plasma is shown. As is seen, our treatment of the in-medium 
strings provides quite reasonable results. The quantity {e — 3p)/T^ turns out to be very 
sensitive to the QCD interactions. Indeed, it is still equal to zero even in the one-gluon- 
exchange approximation provided the temperature dependence of the running coupling 
constant is neglected. The agreement with the lattice calculations could be even better if 
we chose 7 ~ 2/3 like in Ref. P3. Thus, the interesting question arises what additional 



arguments, being able to change the 7- value from 1/3 to 2/3, should be taken into account 
for our picture of the in-medium string screening. In this respect the spatial correlations 
of color charges may be of importance (see Sec. 3). 

In the transchemistry |11] the reduced effective value a = 0.5 GeV/fm = 0.1 GeV^ is 
used which results in Tchem = 185 MeV. It seems that the transchemistry model starts at 
slightly lower temperature where the chemical equilibrium for the quarks would not be 
possible at all. On the other hand, this simulation begins with a huge oversaturation of 
the quark number, so a later reheating of the system brings the massive quark matter in an 
over-critical state. Eventually, expansion and cooling leads to dynamical hadronization 
at low temperatures (T < Tchem), where the quark component cannot be in chemical 
equilibrium any more. 

Now, returning to the Boltzmann statistics, for the critical number density we have 

nchem = (^ a)^/^ . (81) 

9e 
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Substituting this quantity into Eq. 
pressure, 



up and taking ■j = 1/3, A = (2/3)cr, we get a negative 



Pchem 



2X_ 

"sTe' 



-a 



12) 



at the critical point. It means that the mechanical equilibrium ceases at a somewhat 
higher temperature than the chemical one. 

Equations similar to Eq. ([76| ) were considered by Boal, Schachter, Woloshin [|T^ and 
by Moskalenko and Kharzeev [17], as well. However, investigation of the quark plasma at 
zero baryon density, presented in the latter paper, did not take into account the important 
background term $. As to the former one, only the boundary for the high-temperature 
QCD phase, rather than the full thermodynamics, was studied without any reference to 
the problem of thermodynamical consistency. This is why one of the considered spectrum 
of unbound partons in this paper is not consistent with Eqs. (21) and (22). In addition, 
none of these papers uses the relevant classical approximation providing analytical results 
like Eqs. (|79|) , (|8ll) and (^), which would have significantly simplified understanding. 
At last, an advantage of our treatment is that it is based on the elaborated model of 
the in-medium string formation that enables us to derive the mean-field term (|5l|) rather 
than postulate it invoking the nearest-neighbor approximation inspired by the analogy 
with the Ising model (see the third paper in Ref. |ll5|). 



4.5 Fermions with string interaction 

Another interesting case is massless quarks at zero temperature. The number density 
integral is given by 

n = ^ [ e{l2--an-^/^ ~k) k'^dk . (83) 

Here d is the color, light flavor and spin degeneracy factor, Q{x) denotes the Heaviside 
step function. Expression ( p3D can be readily rewritten as 

with C = (d/6TT^y^^. In the situation considered we have one conserved charge: the 
baryon number with the density Uf, = n/3. Then, the chemical equilibrium is specified 
by the relation /i = /i;,/3, where fib is the baryon chemical potential (see Sec. 2). The 
magnitude of fib is determined by the equation 

lOb — ■ 

Using z = (Ufj ^/'^/33/4 



" 81 ~ 3V3 ) ■ (84) 



1 ~ 

ji — ai z (85) 



z 



with (Ti = 6cr/C. Equation (p5D has a solution provided 

l^b > /ichem = 2v^, (86) 
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i.e. the chemical potential (Fermi energy) is larger than the minimum value of (aiz + l/z). 
It means that the Fermi energy of quarks should be larger than l^/oxj?). A typical 
numerical value, /Xchem ~ 2.442 GeV, can be found using d = 12 and a = 0.18 GeV^. 

At finite temperature we obtain a T- and /ife-dependent consistency equation, which 
can be solved only numerically. Fig.^ shows the resulting boundary in the T - fib plane. 

4.6 Chemical ofF-equilibrium in the classical approximation 

If an isolated system is out of chemical equilibrium and expands as a perfect fluid, the 
relation 

dE + p dV = T dS + dNi = (87) 

i 

is fulfilled and the entropy production rate 

i 

is either positive or zero. This means that in the one-component case for quasiparticles 
with positive /i, the corresponding particle number = Vn decreases while it increases 
for negative values of fi. 

Fig.| shows the off-equilibrium chemical potential scaled with the temperature, /i/T, as 
a function of the scaled density ra^/^/T for a one-component, massless Boltzmann gas made 
of particles with the SU (3)-gluon degrees of freedom and interacting via strings. Chemical 
equilibrium corresponds to yU = 0, which is not reachable below a certain temperature. 
Then the strings pull the charges together never reaching a screened equilibrium state: 
the chemical potential remains positive driving the density of this component towards 
zero. 

The situation is more complicated in a many-component mixture due to possible con- 
stituent exchange between different species of quasiparticles. In both cases the chemical 
equilibrium, corresponding to /ij = 0, is stable. In some special cases for non-ideal EoS it 
may happen that the fi{n, T) curve for a constant T (isoterm) does not cross the /i = 
line at all, i.e. no chemical equilibrium is possible and the system is driven towards a 
state with either zero or infinite particle numbers. In a many-component system it means 
that this particular component will dominate or vanish in the mixture. 

Another remark concerns with the chemical potential assigned to the conserved charges 
{e.g. baryon number). This is a physically different situation when the term —fibQb (see 
Sec. 2) is added to the Hamiltonian, which is not compensated by its expectation value in 
the background field. As a consequence, the chemical equilibrium point (if any) is placed 
not at /i = but aX ji = hb- This situation is quite accustom in nuclear physics. 



5 Conclusions 

A useful representation of the conditions of thermodynamical consistency of quasiparticle 
description, Eqs. (^Tj) and (^), has been found. The advantage of this representation 
is that it directly involves the effective quasiparticle spectra, which results in important 
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restrictions to the form of these spectra. In particular, two essential findings can be 
mentioned. If the interaction with surrounding matter is taken into account by introducing 
a mean field, the latter should be either temperature independent (Example 1, Sec. 2) or, 
when in-medium effects are included into the Hamiltonian by means of the effective mass 
(Example 2, Sec. 2), this mass should depend on the temperature or/and the quasiparticle 
density exclusively through the scalar density of quasiparticles. On the large market of 
available phenomenological models these general restrictions in majority of cases were 
used intuitively, but sometimes were erroneously missed. 

The structure of the thermodynamical potential, derived from the medium- dependent 
Hamiltonian in a thermodynamically consistent way, has been further detailized by im- 
plementing the string picture for the interaction between generic constituents. With this 
aim the elaborated mean-field model of the in-medium string interactions has been devel- 
oped. This model supports the use of the inverse power of the color charge density in the 
color mean field (see Eq. (^T])) that was introduced earlier by various authors assuming 
validity of the nearest neighbour approximation [1^ . Results of our treatment of 



the in-medium strings are found to be in reasonable agreement with the lattice data on 
QCD thermodynamics. Further probing of these interactions can be an application of the 
developed equation of state to (hydro) dynamical calculations allowing direct comparison 
with observables. Some steps towards this direction have been done recently by analyzing 
the excitation function for nucleon directed flow fS^ and the relation of the 'softest point' 
of equation of state with chemical freeze-out in heavy ion collisions. 

Along with other results, we would like to comment on the excluded volume modifi- 
cation of the single particle energy. In our treatment the conditions Ct = 0, Ci = have 
no solution in this case without some additional assumptions. A resolution of this issue 
has been done in Ref. and in the first paper of Ref. [O. 
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Figure Captions 



Fig.l. Normalized energy density and pressure of a massless gluon gas witli string-like 
interaction. Dashed line demonstrates the Stephan-Boltzmann regime corresponding to 
the case of a gas of noninteracting gluons. 

Fig. 2. The interaction measure {e — 3p)/T^ of the SU{3) gluon system: circles are 
our data for the massless gluons interacting via screened strings; squares show the lattice 



results from Ref. [32]. Our data are plotted for the case T^ec = 303 MeV. 



Fig. 3. The critical curve on the temperature-chemical potential plane, below which 
chemical equilibrium ceases for a massless quark gas with string-like interaction. 

Fig. 4. The scaled chemical potential jj/T as a function of n^^^ /T for massless Boltz- 
mann gas with string-like interaction. The chemical equilibrium condition is = 0. 
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Figure 1: Normalized energy density and pressure of a massless gluon gas with string- like 
interaction. Dashed line demonstrates the Stephan-Boltzmann regime corresponding to 
the case of a gas of noninter acting gluons. 
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Figure 2: The interaction measure {e — 3p)/T^ of the SU{3) gluon system: circles are 
our data for the massless gluons interacting via screened strings; squares show the lattice 
results from Ref . . Our data are plotted for the case Tjec = 303 MeV. 
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Figure 3: The critical curve on the temperature-chemical potential plane, below which 
chemical equilibrium ceases for a massless quark gas with string-like interaction. 
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Figure 4: The scaled chemical potential fi/T as a function of n}/^ /T for massless Boltz- 
mann gas with string-like interaction. The chemical equilibrium condition is = 0. 
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